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The attractor dimension at the transition to complete synchronization in a network of chaotic units 
with time-delayed couplings is investigated. In particular, we determine the Kaplan- Yorke dimension 
from the spectrum of Lyapunov exponents for iterated maps and for two coupled semiconductor 
lasers. We argue that the Kaplan- Yorke dimension must be discontinuous at the transition and 
compare it to the correlation dimension. For a system of Bernoulli maps we indeed find a jump in 
the correlation dimension. The magnitude of the discontinuity in the Kaplan- Yorke dimension is 
calculated for networks of Bernoulli units as a function of the network size. Furthermore the scaling 
of the Kaplan- Yorke dimension as well as of the Kolmogorov entropy with system size and time 
delay is investigated. 



INTRODUCTION 

Networks of identical nonlinear units which are coupled 
by their dynamic variables can synchronize to a common 
chaotic trajectory [TJ|2]. This phenomenon is of funda- 
mental interest in nonlinear dynamics, with applications 
in neural networks, coupled lasers, electronic networks 
and secure communication (3HS]. For many applications, 
the coupling signals are transmitted with a time delay 
which is much larger than the internal time scales of the 
individual units. A network of identical units with iden- 
tical delay times can synchronize to a common chaotic 
trajectory without any time shift between these units 
(zero lag synchronization) [7J. This can only occur for 
so-called weak chaos [5], i.e., if the largest Lyapunov ex- 
ponent (LE) of the network decays with the inverse delay 
time. 

A chaotic system with time delayed couplings (includ- 
ing self-feedback for a single unit) has special mathemat- 
ical properties. The system becomes infinite dimensional 
and has a continuous spectrum of LEs. The dimension 
of the chaotic attractor increases proportional to the de- 
lay time [HI [ID] • In this contribution we investigate the 
attractor dimension close to the transition to chaos syn- 
chronization. 

The dimension of a chaotic attractor is not uniquely 
defined [IT]. We consider two attractor dimensions: The 
Kaplan- Yorke (KY) dimension [TJ] , Dk y , which is deter- 
mined from the spectrum of the LEs, and the correlation 
dimension, Dc [El [14]. The KY dimension is conjec- 
tured to be identical to the information dimension, Dj. 
This is known as the so-called KY conjecture [T51 IPd] . 
The information dimension in turn is an upper bound 
for the correlation dimension, Di > Dc [17] , 

In this article we study these two attractor dimensions, 
Dky an d Dc, for networks of identical nonlinear units 



coupled by their time-delayed variables. The delay time 
is much larger than the internal time scales, and chaos 
is generated by the coupling. Increasing the coupling 
strength, the system has a transition to complete syn- 
chronization [TS] . 

We show that the KY dimension is discontinuous at 
this transition, it jumps from a high value for the un- 
synchronized to a low value for the synchronized chaotic 
attractor. This is a general results which holds for any 
networks with time-delayed couplings. We numerically 
calculate the KY dimension for networks of iterated maps 
and for two coupled semiconductor lasers. For networks 
of Bernoulli maps we are also able to obtain analytical 
results. As a cross check the correlation dimension of the 
systems is computed as well and (at least) for networks 
of Bernoulli maps it displays a clear discontinuity at the 
synchronization transition in agreement with the jump in 
the KY dimension. Finally we calculate the jump in the 
KY dimension as a function of the size of the network. 

DISCONTINUOUS KAPLAN- YORKE 
DIMENSION 

Let us consider a network of N identical nonlinear 
units. Each unit j = 1, . . . , N has a set of dynamic vari- 
ables Xj (t) which obey the following differential equations 
with time-delayed couplings 

JV 

ij(t) = F[ Xj (t)] + o-Yl G ^ H ^ - T )l ' (!) 

k=l 

The function F describes the local dynamics of each 
unit, the function H couples the time-delayed variables of 
the connected units and the adjacency matrix Gjk defines 
the graph of the network. We restrict this matrix to 
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have only nonnegative entries and a constant row sum 
J^k^jk — 1- Thus the eigenvalue of G with the largest 
modulus is always 71 = 1. The parameter a determines 
the strength of the coupling and r is the delay time of 
the coupling. 

By construction, the synchronization manifold (SM) 
Xj(t) = s(t) is a solution of this network given by the 
equation 



s(t) = F[s(t)} + aH[s(t - t)} 



(2) 



We consider only networks where this equation has 
chaotic solutions for sufficiently large values of er, i.e., 
the dynamics on the SM is chaotic and has at least one 
positive LE. The stability of the SM can be determined 
using the master stability function |19j . An infinitesi- 
mal perturbation of the synchronized trajectory can be 
decomposed into eigenvectors of the coupling matrix G 
with corresponding eigenvalues 7^, k = 1,...,N. The 
amplitudes of the perturbations along these eigen- 
vectors follow the equation 



Ut) = DF[s(t)] + a lk DH [s(t - r)] &(t - 



(3) 



Hence the perturbations are governed by a linear differen- 
tial equation with time-delayed feedback and time depen- 
dent coefficients. For each eigenvalue 7^ of G this equa- 
tion yields a whole spectrum of LEs. A perturbation cor- 
responding to the largest eigenvalue, 71 = 1, describes a 
perturbation within the SM. Hence 71 is called transver- 
sal eigenvalue. In the chaotic regime, the linearized dy- 
namics is unstable for 71 = 1, i.e., it has positive LEs. 
Any arbitrary perturbation has in general components in 
the other eigenmodes with eigenvalues 7fc,fc = 2,...,N. 
These are the so-called longitudinal eigenvalues. Thus 
the SM is stable if eq. ^ yields only negative LEs for all 
transversal eigenvalues such that perturbations transver- 
sal to the SM are decaying exponentially fast. 

We consider only networks where the delay time t is 
much larger than any other time scale of the system [5D] . 
In this case the condition for stability of the SM is given 
by the following equation [5] 



I72I < exp(-A m r) . 



(4) 



Hence the transversal eigenvalue of G with the largest 
absolute value, I72I, and the maximum (longitudinal) LE, 
A,„, of a single unit with feedback, given by eq. de- 
termine the stability of complete chaos synchronization. 
Note that A m depends on the coupling strength a, there- 
fore eq. Q determines the critical coupling strength a c 
where chaos synchronization appears. 

In this work, we change a control parameter of the 
system such that the system exhibits a transition from 
from a synchronized to an unsynchronized state. At this 



transition we investigate the change in the attractor di- 
mension. 

A quantitative measure for the structure of the attrac- 
tor is the KY dimension. It is defined by the spectrum of 
LEs obtained from eq. Considering a discrete spec- 
trum Ai > A2 > ■ ■ ., the KY dimension is defined as the 
largest number M for which the sum of LEs is still pos- 
itive plus an interpolation term which yields the fractal 
part of the dimension, 



D KY = M 



M+l| 



(5) 



Now we want to show that the attractor dimension at 
the transition to chaos synchronization is discontinuous. 
The LEs are given by eq. ^ where one obtains a whole 
spectrum of LE for each eigenvalue 7&. 

First, consider stable synchronization. For the lon- 
gitudinal eigenvalue 71 = 1 the perturbations have in 
general positive as well as negative LEs. For all other 
(i.e., transversal) eigenvalues the corresponding LEs are 
negative since the SM is stable. But close to the synchro- 
nization transition the maximum LE of (at least) one of 
these transversal Lyapunov spectra approaches the value 
zero, as shown in Fig.JT] Consequently, close to the tran- 
sition these spectra contribute to the definition of the KY 
dimension, defined by eq. ([5]). 

However, this cannot be true. In case of stable syn- 
chronization the trajectory of the network, determined 
by eq. pi), is completely restricted to the SM. Two 
neighboring trajectories inside the SM deviate from each 
other according to eq. ([3| with 71 = 1. Any transver- 
sal eigenvalues cannot contribute to the dynamics inside 
the SM. Thus, only the longitudinal LE spectrum can 
contribute to the KY dimension, i.e., the sum in eq. §5§ 
must only run over the LEs of the 71 spectrum. 

Now, consider the dynamics outside but still close to 
the synchronization transition. In this case, the pertur- 
bations cannot be decomposed into the eigenmodes of the 
coupling matrix G since the coefficients of the linearized 
equations of eq. depend on each single node. But 
close to a supercritical transition we expect that the LEs 
are continuous as a function of a and the structure of the 
LEs is still similar to the one inside the SM, as it is the 
case in Fig. [T] 

In the desynchronized regime all eigenvalues contribute 
to the KY dimension. Consequently, the KY dimension 
must be discontinuous at the transition to chaos synchro- 
nization. 

Chaos synchronization can only occur for weak chaos, 
where the largest LE scales with 1/r. In this case the 
spectrum of LEs is dense. Hence, a large number of ad- 
ditional LEs contributes to the KY dimension at the tran- 
sition, and we expect the jump of the KY dimension to 
be of the order of the size, N, of the network. 
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Figure 1. (Color online) Lyapunov spectra vs. coupling 
strength e for a system of three all-to-all coupled Bernoulli 
maps with the parameters a = 1.5 and r = 20. The blue 
crosses show the longitudinal spectrum associated with 71 = 1 
and the red dots show the transversal spectrum associated 
with 72 = 73 = —1/2. Both spectra were obtained by solving 
the polynomial equations derived from the master stability 
function. The black solid lines show the Lyapunov spectrum 
computed by Gram-Schmidt from simulations of the full sys- 
tem. The vertical dashed line indicates e c where to its right 
the system is synchronized. 

From the Lyapunov spectrum we can also calculate the 
Kolmogorov entropy which quantifies the predictability 
of the system [5T] . It is defined as the sum over all posi- 
tive LEs 

A'^^A,, forA;>0. (6) 

i 

Since only positive LEs contribute, the entropy is always 
defined with the complete spectrum of LEs, hence it does 
not show a jump at the transition to chaos synchroniza- 
tion. Nevertheless, at the transition at least one band of 
LEs adds to the entropy and we expect a discontinuous 
derivative of K (a) 

The jump in the dimension of the chaotic attractor and 
the kink in the entropy are general results which should 
hold for any chaotic network at the transition to chaos 
synchronization. In the following section we calculate the 
attractor dimensions T)ky & n d Dq and the entropy K 
for networks of iterated maps. 



ITERATED MAPS 

The previous general statement holds not only for dif- 
ferential equations but also for networks of iterated maps 
with time-delayed coupling. Since such models are easier 
to analyze than continuous systems we investigate net- 
works of iterated maps in this section. Each unit j has 
a one-dimensional variable in the unit interval x\ S [0, 1] 



which is updated in discrete time steps t according to the 
following equation 

4 + x = (1 - e)f(4) + ^G ]k f{x k t _r) ■ (?) 

k 

The parameter e is the coupling strength, but since 
the dynamic variable should stay in the unit interval, we 
subtract the undelayed term with strength e. For the 
function f(x) we use the Bernoulli shift and the asym- 
metric tent map, 

Bernoulli map 

x t+ i = (ax t ) mod 1 

Tent map 

( i x t for < x t < a 

Xt+1 = \ a ^- a {l-x t ) fora<x t <l 

The Bernoulli map is chaotic for a parameter a > 1. 
For the chaotic tent map the parameter a is chosen such 
that the value Xt+i stays in the range [0, 1], i.e., < a < 
1. For maps, we have a discrete delay time r. As before, 
the coupling matrix G has constant row sum ^2 k Gjk = 1. 

The synchronized state is a solution of these equations 
and reads as 

a t+1 = (1 - e)f(s t ) + e/(a t _ T ) . (8) 

As before, the perturbations of the SM can be associ- 
ated with the eigenvalues 7^ of the coupling matrix G. 
The amplitude £ t fe of each mode obey the linear master 
stability function determined by 

= (1 - e)f(s t ) +e7k f(s t . T ) g_ T . (9) 

Note that for the Bernoulli network the derivative /' = 
a is constant, therefore one only has to analyze linear 
equations with constant coefficients. For the tent map, 
however, the derivative can take on two different values 
and, hence, the coefficients change with time. 

Since r is discrete, for each mode with eigenvalue 74, 
one obtains r + 1 many LEs. Fig.[l]shows an example for 
a triangle of all-to-all coupled Bernoulli units with 71 = 1 
and 72 = 73 = —1/2. For all values of e the system is 
chaotic, since the largest LE of the 71 band is always 
positive. The transition to chaos synchronization occurs 
at e c where the maximum LE of the 72,3 band crosses the 
value zero. 

The Lyapunov spectrum is in general obtained from 
a Gram-Schmidt orthonormalization procedure accord- 
ing to Farmer [TU]. The system's equation ([!]) and ([7|, 
respectively, are linearized around the chaotic trajectory 
and simulated for a set of orthogonal perturbation vectors 
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which have to be re-orthogonalized after an appropriate 
amount of time. The Lyapunov spectrum is computed 
from the change in magnitude of the perturbation vec- 
tors. 

For Bernoulli networks we can also derive a polynomial 
equation of degree r + 1 for the different eigenvalues jk of 
the adjacency matrix G from which the Lyapunov spec- 
trum can easily be calculated. The polynomial equation 
reads as follows 



(1 - e) az r + eaj k 



(10) 



where the LEs are given by A = In | z\ [22] . This equa- 
tion still holds in the desynchronized region since it does 
not depend on the trajectory of the system. Thus the LE 
spectrum can be calculated from eq. (|9| and ( 10 ), respec- 
tively, in the complete parameter space. Note that both 
methods - orthonormalization procedure and polynomial 
equation - which compute the spectrum in completely 
different ways yield the same results, see Fig. [T] 

From Figs. [T] and [3] it can be seen that for e- values close 
to 1 the LEs cluster into bands. This can be understood 
as follows. For e = 1, the dynamical equations Q are 
given by 



N 

fc=i 



G jk f(xt T ) 



(11) 



Since the state at time t + 1 is only influenced by the 
state at time t — r, the system is effectively given by 
r + 1 uncoupled identical systems of the form 



E G ^/(^ fe ) 



(12) 



Each of these r + 1 systems is iV-dimensional and gives 
rise to N Lyapunov exponents. Since we have r + 1 iden- 
tical systems, each of these N exponents is r + 1 times 
degenerate. This holds as long as each effective system 
evolves on the same chaotic attractor, and thus does not 
rely on synchronization. 

For e < 1 the first term in eq. Q leads to a coupling 
between these effective systems and thus removes the de- 
generacy. 

It is insightful to discuss this lifting of degeneracy for 
the case of the Bernoulli maps. For e = 1 the solutions of 
the variational eq. ( 10 ) are given by the complex (r + 1)- 
th roots of 7fca 



,(o) _ 



(7feo)-+i e 1 ^ 1 (/ = 0, ...,r) 



(13) 



and the corresponding (r + 1) Lyapunov exponents are 
all equal and are given by 



A (0) =ln|z (0) | =ln| 7fc a|^r 



(14) 



Here the superscript (0) indicates the zero-th order in an 
expansion in 1 — e. 
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Figure 2. (Color online) Perturbation expansion of eq. ( 10 I 
up to first order in /i = (1 — e) (red lines) and exact location of 
eigenvalues (blue dots). Parameters: a = 0.4, 74. = 1, r = 20. 



For e ^ 1, we can make a perturbation expansion of 
eq. ( 10 ) in the small parameter /1 = 1 — e. Inserting the 
ansatz z = + fift 1 ', where (3^ is a coefficient that 



needs to be determined, into eq. ( 10 ) yields the solutions 
up to first order in in \x 



,(0) 



,(0). 



7fe 



1 



(15) 



The values of the corresponding Lyapunov exponents 
(A^ 1 ) = ln|zW|) are depicted in Fig. ^1 as a function of e. 
One clearly sees the lifting of degeneracy due to e < 1. 
Although for general maps / one cannot write down an 



equation such as ( 10 1 in the unsynchronized regime, for 
e = 1 the degeneracy follows rigorously from the discus- 
sion above and for e < 1 the interaction will generically 
lead to a lifting of degeneracy similar to that shown in 
Fig. [2] 

In the synchronized regime the linear equations ^ and 
([9]) can be used to compute the LE spectrum for any it- 
erated map f(x). Hence the LEs can be clustered into 
bands according to the eigenvalues of the coupling matrix 
G. However, for the unsynchronized system these linear 
equations only hold for Bernoulli networks where the co- 
efficients are constant and in particular independent of 
the systems trajectory. In general, in the desynchronized 
regime we need to evaluate the linearized equations of the 
full system, eq. and Q, to obtain the LE spectrum, 
and cannot restrict to the master stability function, eq. 
([3]) and ^9|. But close to the transition we expect the 
spectra obtained from the master stability function to 
approximate the true spectra very well. Surprisingly, for 
the tent map, the results coincide very well not only close 
to the synchronizations transition but for all values of e 
when for the desynchronized system in the master sta- 
bility function, eq. (|9|), the dynamics of a single unit is 
inserted, e.g., s t and s t _ T is replaced by x\ and x\_ r , re- 
spectively. A comparison of the different spectra is shown 
in Fig. [3] In contrast to to Fig. [I] the blue and red lines 
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Figure 3. (Color online) Lyapunov spectra vs. coupling 
strength e for a system of three all-to-all coupled Tent maps 
with the parameters a = 0.4 and r = 10. The blue crosses 
show the longitudinal spectrum associated with 71 = 1 and 
the red dots show the transversal spectrum associated with 
72 = 73 = —1/2 computed by Gram-Schmidt from simula- 
tions of the master stability function. The black solid lines 
show the Lyapunov spectrum computed by Gram-Schmidt 
from simulations of the full system. The vertical dashed line 
indicates e c where to its right the system is synchronized. 
Note that for a better visibility a system with r = 10 is plot- 
ted. 



in Fig. [3] are obtained using the Gram-Schmidt proce- 
dure from simulating the master stability function which 
is strictly only valid for the synchronized regime. The 
black line is obtained using the Gram-Schmidt procedure 
on the full system's equations and therefor yields the cor- 
rect results not only for the synchronized but also for the 
unsynchronized regime. Within the synchronized regime 
the results match up to numerical accuracy, whereas out- 
side of synchronization the results of the two methods 
deviate since the master equation is no longer valid. 

The KY dimension for the triangle of Bernoulli and 
tent maps, respectively, is shown in Fig. [4] The upper 
curve shows the KY dimension when the complete set of 
LEs is used in eq. while the lower curve uses only 
the LEs of the SM which are obtained from simulating 
a single unit. Note that an upper bound of the dimen- 
sion is 3(r + 1), i.e., the full system's dimension, for the 
desynchronized triangle and r + 1 for the manifold. 

As discussed in the previous section, in the synchro- 
nized region the lower curve is valid whereas in the desyn- 
chronized region the upper curve is valid. Thus the KY 
dimension jumps to a lower value when the parameter e 
is increased above e c . 

In the following we consider a pair of maps with self- 
feedback. The dynamic equations of the system read 



(1 - e)f{x\) + enf{x\_ T ) + e(l - k)/(zJ_ t ) , (16) 



50 



40 



30 



20 



10 



(a) 









0.2 



0.4 0.6 

£ 



50 



40 



Q 30 



20 



10 



- (b) 









0.2 



0.4 0.6 

£ 



with i,j € {1,2}. The parameter e is, as before, the 
coupling strength of the delayed terms to the internal dy- 



Figure 4. (Color online) KY dimension Dky with respect to 
the coupling strength e for a system of three all-to-all coupled 
(a) Bernoulli and (b) tent maps, respectively. The upper blue 
line shows Dky of the full system, the lower red line shows 
Dky of the SM. The vertical dashed line indicates e c where 
to its right the system is synchronized. The parameters are 
t — 20, a — 1.5 for the Bernoulli map and a — 0.4 for the 
tent map, respectively. 



namics and the parameter k determines the ratio between 
the external coupling and the self-feedback. The synchro- 
nized trajectory of the system also follows eq. ^ which 
does not contain the strength n of the self-feedback. 
Hence the synchronized trajectory, i.e., the SM is inde- 
pendent of k and only changes with the coupling strength 
e. Fig. [6] shows the KY dimension as a function of k. At 
the transition the dimension jumps from the upper (blue) 
curve to the lower (red) constant value. 

According to the discussion of the previous section the 
KY dimension has to jump at the synchronization tran- 
sition. This qualitative prediction of the KY conjecture 
should be valid for any measure of the dimension of the 
chaotic attractor. Thus, we also computed the correla- 
tion dimension of the system to compare it to the KY 
dimension and to check whether the dimension indeed 
jumps as we argue. For this purpose we analyzed the 
system's trajectories, i.e., the time series of the system 
using the TISEAN package of Kantz and Schreiber [33] . 
In particular, the correlation function C(£) according to 
Grassberger & Procaccia was computed which scales as 
a power law C(£) oc £ Dc with the exponent being the 
correlation dimension [T2J [H] . A straight line was fitted 
to different correlation functions of different embedding 
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Figure 5. (a) Correlation function C(£) and (b) local slope 
for different embedding dimensions m (different curves) 
computed from a timeseries of a system of two mutually cou- 
pled tent maps. The time delay is r = 5, the length of the 
time series is I = 10 6 and the other parameters are a = 0.4, 
e — 0.45 and re = 0.35. Vertical dashed lines indicate the 
range which was used for the fit. 



dimensions in a double-logarithmic plot and at the same 
time the results were cross-checked in plots of the local 
slopes of the correlation function d(£) = <9C(£)/<9£ in 
which the power law behavior corresponds to a plateau. 
For more details on how to actually compute the corre- 
lation dimension the reader is referred to [22 IMj ■ Note 
that this method allows a reliable calculation of the cor- 
relation dimension for small values of the delay t, only. 

A typical plot of the correlation functions as well as 
the local slopes is shown in Fig. [5] This figure shows 
that the extrapolation of the slope to low values of the 
radius £ is difficult. Hence our values for the correlation 
dimension give a lower bound to the actual correlation 
dimension since due to the limited computational power 
it was not possible to analyze very small values of the 
distance £ with an appropriate accuracy. 

Fig. [6] shows the results for the KY and correlation di- 
mension as a function of re. The KY dimension is larger 
than the correlation dimension in agreement with known 
theoretical inequalities. For the Bernoulli map the corre- 
lation dimension displays a clear jump at the transition to 
synchronization, in agreement with our theoretical pre- 
diction. For the tent map, however, the discontinuity 
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Figure 6. (Color online) Attractor dimension vs. coupling 
strength re for a system of two mutually coupled (a) Bernoulli 
and (b) tent maps, respectively. The time delay is r = 5. 
The upper blue line shows Dky of the full system, the lower 
red line shows Dky of the SM corresponding to 71 and the 
green curve in between shows Dc- The vertical dashed lines 
indicate re c . The parameters are r = 5, e = 0.45, a — 1.5 for 
the Bernoulli map and a = 0.4 for the tent map, respectively. 

is not clearly visible from our results. As stated above, 
for small distances £ the correlation function C(£) shows 
large fluctuations due to the limited statistics. From our 
results of Fig. [5] we cannot rule out that the local slopes 
might not be saturated yet. Thus the obtained results 
only serve as a lower bound which, according to our re- 
sults, seems to increase with longer trajectories and bet- 
ter statistics. Consequently, the results of Fig. [6] do not 
rule out a discontinuous behavior of the attractor dimen- 
sion. In any case, the synchronization transition is clearly 
visible in the discontinuous slope of the correlation di- 
mension. 

The Kolmogorov entropy, computed from eq. ([6]), is 
shown for a pair of Bernoulli and tent maps, respectively, 
in Fig. [7j In the synchronized region only the 71 band 
has positive LEs which contribute to the Kolmogorov en- 
tropy. At the synchronization transition a kink in the 
entropy as a function of feedback strength can be seen 
when suddenly LEs from the other band contribute. 

The result for the attractor dimensions and the pre- 
diction entropy are very similar for Bernoulli and tent 
maps. But we found a qualitative difference between the 
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Figure 7. Kolmogorov entropy K vs. coupling strength « for 
a system of two mutually coupled (a) Bernoulli and (b) tent 
maps, respectively. The vertical dashed lines indicate k c . The 
parameters are r = 20, e = 0.6, a — 1.5 for the Bernoulli map 
and a = 0.4 for the tent map, respectively. 



two models for the cross-correlations, C, and the syn- 
chronization probability, </>, of a system of two mutually 
coupled units. The synchronization probability measures 
the fraction of time where the two trajectories are closer 
than some threshold O [25] . Fig. [8] shows C and <f> with 
respect to the coupling strength k for fixed e. At the 
critical coupling k c , i.e., at the synchronization transi- 
tion both quantities, C and <f>, jump from a very low 
level to complete synchronization, C = <j> = 1, for a sys- 
tem of Bernoulli maps, whereas for the tent map C and 
(j> increase continuously to C = <j> = 1. 

The numerical results indicate that the synchroniza- 
tion transition for the tent map is of a supercritical type. 
That is, close to the transition to synchronization there 
is a stable trajectory close to the SM and the dynamics 
is nearly synchronized. We thus observe a smooth tran- 
sition to synchronization. In the Bernoulli map, on the 
other hand, the transition is of a subcritical type. Note 
that in both cases, a jump of KY dimension is predicted 
due to the bands of LEs. 



COUPLED LASERS 

An important application of eq. is the modeling 
of semiconductor lasers which are coupled by their mu- 
tual laser beams. To a good approximation, the dynam- 
ics of the laser intensity can be described by the Lang- 
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Figure 8. (Color online) Cross-correlation C (blue squares) 
and synchronization probability (j> (red dots) vs. coupling 
strength k for a system of two mutually coupled (a) Bernoulli 
and (b) tent maps, respectively. The step size is Ak = 10~ 3 
and the other parameters are r = 20, e = 0.6, a = 1.5 for the 
Bernoulli map and a = 0.4 for the tent map, respectively. The 
threshold for <f> was set to = 0.01, i.e., the trajectories were 
assumed to be synchronized when they were closer together 
than 1% of their maximum distance. 



Kobayashi rate equations [26 . The Lang-Kobayashi 
equations describe the dynamics of a laser with delayed 
feedback (or delayed coupling) in terms of a slowly vary- 
ing complex electric field E(t) and a population inversion 
n(t). For our network of coupled lasers, the correspond- 
ing equations of Lang-Kobayashi type in dimensionless 
form are given by 

E\i) = hi + ia^i^E^t) + aJ^GtjE 1 (t - t) 



Til 1 (t)=p- n l (t) - [1 + n l (*)] | E l (t) \ 



(17) 



where E l {t) is the envelope of the complex electric field 
and n l (t) is the renormalized population inversion of the 
charge carriers of laser i. The model parameters are 
summarized in Table [Tj The dimensionless delay time of 
t = 100 translates to a delay time of the order of magni- 
tude Ins. 

A network of coupled lasers modeled by the Lang- 
Kobayashi equations can be written in the form of eq. ([!]) , 
where x l (t) = (n l , Re E l , Im E l ) is now three-dimensional 
and contains the real and imaginary part of the electric 



Parameter 



Symbol Value 



Linewidth enhancement factor a 4 
Time-scale separation of 

carrier and photon lifetimes T 200 

Injection current p 0.1 

Coupling strength a 0.12 

Coupling delay time r 100 

Table I. Parameters for the simulation of the Lang-Kobayashi 
equations. 



field E l and the charge carrier inversion n l of the i-th 
laser. A single laser is not chaotic, but the delayed feed- 
back and/or coupling renders the system chaotic. The 
linear coupling function H is represented by the matrix 
H = (oioj, corresponding to all-optical coupling as in 



Eq. (17) 



We consider a pair of lasers with overall coupling 
strength a and self-feedback strength k. That is, the 
coupling matrix is given by 



G = 



K 1 — K 
1 — K K 



(18) 



Note that in the SM, the trajectory eq. |2]) does not de- 
pend on the parameter n. Similar to the case of the maps 
discussed before, the spectrum of LEs is obtained from 
a Gram-Schmidt orthonormalization procedure accord- 
ing to Farmer [10 . From this spectrum we obtain the 
KY dimension using eq. ^. Fig.^a) shows D KY as a 
function of k. The dashed blue curve was obtained using 
the complete spectrum, while the solid green curve uses 
only the spectrum inside the SM. The vertical gray lines 
denote the boundary of stable synchronization of the two 
lasers with the given parameters; synchronization is sta- 
ble between the two lines. 

Outside the synchronization region the complete spec- 
trum has to be used when computing the KY dimension. 
However for chaos synchronization only the spectrum in 
the SM, eq. needs to be used. Thus, at the transition 
the KY dimension has to jump from a high (dashed blue) 
to a lower value (solid green) for the SM. 

Unfortunately, we are not able to calculate the correla- 
tion dimension of the laser rate equations. Due to the de- 
lay term which makes the system high-dimensional, the 
available algorithms do not produce reliable results, to 
our knowledge. 

We have also calculated the Kolmogorov-Sinai entropy 
K. It is denned as the sum of all positive LEs. Fig. [9jb) 
shows the result K(k) for the pair of lasers, as before. 
Again, the dashed blue curve was obtained using the com- 
plete spectrum, while the solid green curve uses only the 
spectrum inside the SM. At the transition, the derivative 
of K(k) is discontinuous. In the desynchronized region, 
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Figure 9. (Color online) (a) Kaplan- Yorke dimension Dky 
and (b) Kolmogorov-Sinai entropy K of two coupled semi- 
conductor lasers in dependence on the relative self-feedback 



strength k, cf. Eq. (181. The dashed blue and solid green 
lines are obtained using the complete spectrum and the using 
the spectrum inside the SM only, respectively. 



the second band of LEs crosses zero, as shown in Fig. [T] 
and|3]for the Bernoulli and the tent map triangle, respec- 
tively, and therefore it contributes to the Kolmogorov- 
Sinai entropy. 



NETWORKS 

In this section, we investigate the transition to chaos 
synchronization in large networks. For simplicity, we 
consider Bernoulli networks of N units with all-to-all 
couplings without self-feedback. The coupling matrix 
G of eq. |7| has the eigenvalues 71 = 1 and jj = 
-l/(N-l),Kj<N. 

Fig. [TO] shows the spectrum of LEs as a function of e for 
a system of five all-to-all coupled Bernoulli units. Since 
the eigenvalue corresponding to the transversal spectrum 
has the multiplicity four this spectrum is four-fold degen- 
erated. For e — 0, the uncoupled Bernoulli units have the 
LE A = log a. Hence the KY dimension is Dky = N, i.e., 
the full phase space. As discussed before, at the transi- 
tion to synchronization only the 71 band contributes to 
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Figure 10. (Color online) Lyapunov spectra vs. coupling 
strength e for a system of five all-to-all coupled Bernoulli 
maps with the parameters a = 1.5 and r — 20. The blue 
dashed lines show the longitudinal spectrum associated with 

71 = 1 and the red solid lines show the transversal spec- 
trum associated with the four-fold degenerated eigenvalue 

72 = —1/4. The vertical dotted line indicates the transition 
between strong and weak chaos and the vertical dashed line 
indicates e c where to its right the system is synchronized. 



Dky, hence the KY dimension jumps at the transition. 
The entropy has a kink, since at the transition the N — 1 
many 7j bands do not contribute any more to the en- 
tropy. 

We consider the limit of large delay times, r — > oo. 
In this limit, using eq. 



(10) 



we find the critical value at 
which the transition to synchronization occurs to be e c = 
(a — 1)/ (<z(l — |7j|))- For a coupling strength e < 1 — 1/a 
the system is in the regime of strong chaos otherwise in 
the regime of weak chaos [5]. In general the transition 
from strong to weak chaos does not coincide with the 
synchronization transition which usually occurs at larger 
values e c > 1 — 1/a, see eq. Q. Only in the limit of 
N —> oo both transitions fall together. The transition 
from strong to weak chaos is defined by the change in 
sign of the so-called local Lyapunov exponent [8]. For 
strong chaos it is positive and for weak chaos negative. 

For strong chaos the maximum LE is of order one 
whereas for weak chaos it scales as 1/r. Since each mode 
of the network has r many LEs of the order of 1/r, the 
KY dimension increases linearly with the delay time r 
whereas the Kolmogorov entropy K is independent of r. 
For a large enough coupling strength, e, such that the sys- 
tem is synchronized, the dimension is determined solely 
by the SM and therefore the dimension is independent of 
the number of units, N, and cannot exceed a value larger 
than the delay time of a single unit. 

Fig. [ll] shows the KY dimension and the Kolmogorov 
entropy as a function of the system size, N, in the regime 
of strong chaos for different delay times r. It is clearly 
visible that for large delay times the different plots of 
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Figure 11. (Color online) (a) KY dimension Dky and (b) 
Kolmogorov entropy K with respect to the number of all-to- 
all coupled Bernoulli maps TV for different delay times t. It 
starts from r = 80 (uppermost curve) and increases in steps 
of At = 40 up to t = 200 (lowermost curve). The other 
parameters are a — 1.5 and e = 0.2. 



Dky/tN corresponding to different delays nearly co- 
incide. Thus the KY dimension scales linearly with r 
whereas the Kolmogorov entropy is independent of r. 
Both, entropy as well as dimension increase with sys- 
tem size N and it seems that for large N both quantities 
scale linearly with N. 

As discussed before, the transition to synchronization 
occurs at e c which depends on the number of units in 
the network. In the limit of large delay times r — > oo 
we find e c = (a — l)/(a(l — |7j|))- Thus e c decreases 
monotonically towards e = 1 — 1/a for TV— > oo. At 
the transition only the 71 band contributes to Dky and 
hence the KY dimension jumps. 

Fig. [12] shows the jump of the KY dimension at the 
transition as well as the critical coupling strength e c at 
which the jump occurs as a function of network size N for 
different delay times. Since e c depends on N, the jump 
has a nonmonotonic behavior. For large enough network 
sizes the jump AD scales linearly with the number of 
units. The slope of this linear relation is approximately 
two for small N, AD/N ps 2, and approximately one 
for large N, AD/N ps 1. The transition between the 
two different slopes is related to the dip in the critical 
coupling strength e c (N) which is due to the fact that the 
system is not yet in the limit r — > 00. 
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The slope AD/N depends on the order of limits. If we 
take the limit r — > oo first the jump AD scales linearly 
with the number of units N with a slope of two. If we 
however take the limit N — > oo first it scales with a slope 
of one for any value of r. 

In Fig. [12] (c) the same data as in panel (b) is plotted as 
AD/N versus N/t leading to data-collapse. This shows 
that the jump in the Kaplan- Yorke dimension is in the 
limit of large r determined by the scaling law 

AD NtP(N/t) , 

where ip is the scaling function depicted in panel (c). 

SUMMARY AND DISCUSSION 

Networks of identical nonlinear units can synchronize 
to a common chaotic trajectory. Although the time-delay 
of the coupling can be very large, the units can synchro- 
nize without any time shift. We investigated the transi- 
tion to chaos synchronization in the limit of very large 
delay times. General arguments about the dynamic of 
such a network predict a jump of the KY dimension of 
the chaotic attractor when the network synchronizes. In 
addition, the Kolmogorov prediction entropy should show 
a discontinuous slope. 

We tested these general predictions for networks of it- 
erated maps. For Bernoulli maps, our numerical results 
show a clear discontinuous behavior of the attractor di- 
mension. The KY as well as the correlation dimension 
jump to a low value when the network synchronizes. For 
tent maps, the numerical results of the correlation dimen- 
sion are not so clear due to large statistical fluctuations 
caused by limited computational power. Nevertheless, 
our results indicate a jump in the attractor dimension, 
too. In both cases the prediction entropy shows a dis- 
continuous slope. 

The KY dimension was also calculated for the rate 
equation of semiconductor lasers. Again, our general ar- 
guments give a discontinuous KY dimension at the tran- 
sition. Unfortunately, we were not able to calculate the 
correlation dimension for this case. 

For Bernoulli networks we numerically calculated the 
KY dimension as a function of system size N and delay 
time r in the region of strong chaos. The dimension scales 
with Nr. The jump of the dimension at the transition 
scales with N, as well. 

Our results show that one has to use the standard def- 
inition of the KY dimension, eq. ([5]), with care. If we use 
eq. ^ with all possible LEs we obtain a wrong result. 
Here we argued that not all LEs contribute to eq. ([5| in 
case of synchronization. This argument predicts a jump 
in the attractor dimension, in agreement with our nu- 
merical results of the correlation dimension. However, in 
general it may not be obvious which LEs contribute to 
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Figure 12. (Color online) (a) Critical coupling strength e c 
and (b,c) jump of the KY dimension ADky with respect to 
the system size N of all-to-all coupled Bernoulli maps. The 
different curves correspond to different delay times, r. Lowest 
curve is for r = 40 and increases by steps of At = 20 up to 
r = 160 (upper curve). For e c the theoretical value in the 
limit of t — ^ 00 is plotted as well (upper dashed line). The 
parameter of the Bernoulli map is a = 1.5. 



the KY dimension. For example, if we distort the syn- 
chronization manifold by a nonlinear transformation of 
one unit the dimension does not change but we do not 
know which LE have to be omitted in this case. Also, 
in the case of generalized synchronization the dynamics 
is restricted to a low dimensional manifold which rules 
out the majority of LEs. But again, we do not know in 
advance which LEs have to be omitted from eq. ([5]). 

Our argument of omitting bands of negative LEs relies 
on the fact that the dynamics is restricted to the SM. 
However, a tiny detuning of the nonlinear units leads to 
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imperfect synchronization and this argument is no longer 
valid. Thus, the attractor dimension should immediately 
jump to a high value for any amount of detuning. We 
have tried to calculate the dimension in the limit of zero 
detuning but our numerical results did not allow a con- 
clusive statement. 
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